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We propose a mechanism to pin skyrmions in chiral magnets by introducing local maximum of 
magnetic exchange strength, which can be realized in chiral magnetic thin films by engineering the 
local density of itinerate electrons. Thus we find a way to artificially control the position of a single 
skyrmion in chiral magnetic thin films. The stationary properties and the dynamical pinning and 
depinning processes of an isolated skyrmion around a pinning center are studied. We do a series of 
simulations to show that the critical current to depin a skyrmion has linearly dependence on the 
pinning strength. We also estimate the critical current to have order of magnitude 10^ ~ 10*A/m^. 
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I. INTRODUCTION 

Vortices have long been the standard topological ob- 
jects in condensed matter physics. They are topolog- 
ically stable to perturbations, localized in space, and 
they can be moved, created or annihilated. In recent 
years, another kind of topological object "skyrmion" is 
observed in chiral magnets. A skyrmion is a configuration 
of three dimensional unit vectors in a two dimensional 
space which wraps around the unit sphere once when we 
move from the center of the space to infinity. Skyrmions 
were first proposed in high energy physics [1: and consid- 
ered as an important class of excitations in ferromagnets 
in condesed matter physics. Recently they are theoret- 
ically predicted |1] and experimentally observed in the 
ground state configuration of chiral magnets in a mag- 
netic field |3H6]. The skyrmions either appear isolated 
or form a close-packed lattice. The skyrmion lattice is 
also observed in a single layer of metal film [T with each 
skyrmion having atomic size. Furthermore, a multiferroic 
skyrmion lattice in which spin couples to electric dipole 
is reported very recently |H] . The phase diagram and the 
spinwave spectrum of skyrmion lattice in chiral magnets 
are very well understood in the framework of a classical 
Heisenberg spin model with Dzyaloshinsky-Moriya (DM) 
interaction p,[5HT5]. 

A skyrmion carries a topological charge, and in a con- 
ducting chiral magnet this topological charge is a quan- 
tized emergent magnetic fiux in the view of an adiabati- 
cally moving electron. In real materials, the diameter of 
a skyrmion is around tens of nanometers and the aver- 
aged magnetic field is around several Teslas so the 
coupling between the skyrmions and the conducting elec- 
trons is quite strong. It is observed that this emergent 
magnetic field gives rise to anomalous Hall effect signal in 
the A-phase of chiral magnet MnSi |17l ITS] . The counter- 
active of the conducting electrons to the skyrmion is the 
spin-transfer-torque |19ff25] . It is predicted and observed 
that this torque will push the skyrmions to move in the 
same direction of the current [16, 19, 26 . It is shown 
[26] that the kinetic term in the Lagrangian of a single 
skyrmion makes the skyrmion attain a transverse mo- 



tion with respect to the applied forces, and the averaged 
current-skyrmion interaction gives rise to a skyrmion mo- 
tion parallel to the applied electric current. A moving 
skyrmion can be thought of as a moving magnetic fiux 
and a static magnetic field in a moving frame induces an 
electric field in the laboratory frame according to Fara- 
day's Law. This emergent electric field is measured by 
Schulz et. al. very recently |16| . 

Ultimately we want to manipulate these skyrmions in 
artificial devices. Very recenlty, Tchoe and Han showed 
that a circular current can generate single skyrmions 
without the aid of magnetic field Given that a 

skyrmion is created, another important point of manipu- 
lation is to put a skyrmion at the precise place we want it 
to be. So the pinning and depinning effects of skyrmions 
in chiral magnets is an important issue. In this paper, we 
propose a pinning mechanism of skyrmions by introduc- 
ing local maximum of the magnetic exchange strength. 
We study the effect of such a pinning center by both sta- 
tionary analysis by the energy functional and dynamical 
analysis by the Landau-Lifshitz-Gilbert (LLC) equation. 
Then we find a way to manipulate the skyrmions by elec- 
tric current in an artificial pinning center lattice, which 
may be realized by putting patterned metal grains on the 
surface of a chiral magnetic thin film. We orgnize the 
paper as follows. In section II, we give the field theory 
description of chiral magnets and the stationary analy- 
sis. In section III, we formulate the dynamics of chiral 
magnets by LLG equation on a square lattice and show 
the results of the simulations of the dynamical pinning 
and depinning processes. In section IV, we give some 
discussion and further outlook. 



II. STATIONARY ANALYSIS 
A. Chiral magnet with a pinning center 

The energy density for the chiral magnet is |12) : 
^iot^^-^Y.^d,nf + D{r)n-{V^n)-B-n, (1) 
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where n (r) represents the orientation of local magnetic 
moment with n? = 1, /i = x,y,z, J (r) refers to the 
local ferromagnetic exchange strength, D (r) the local 
strength of DM interaction and B the uniform applied 
magnetic field. In previous works by mean-field theories 
and numerical simulations we already know, for a certain 
range of B, that a skyrmion solution is stable pi [T^ [T5] 
and the trial function for a single skyrmion is: 



n (p, 0, z) — sin [6 (p)] (f) + cos [9 (p)] z, 



(2) 



where we have adopted cylindrical coordinate system 
with local unit vectors p, (j) and z. 9 is the angle be- 
tween n and z. 

To study the pinning effect caused by inhomogeneity, 
we let J and D be space dependent. For simplicity, we 
consider an axisymmetric case such that J and D are of 
p-dependence merely. We choose 



J (p) = J (oo) < 1 -I- A exp 



(3) 



to represent a local maximum of magnetic exchange 
strength. The local magnetic exchange strength is de- 
termined by the local density of itinerate electrons |26| . 
Here we let D (p) and J {p) have constant ratio k = 
D/{2J) because the density of itinerate electrons will 
not change the ratio of spin dependent hopping integral 
to spin independent hopping integral in the underlying 
Hubbard model. Introducing /3 = B/[2J{oo)], we can 
rewrite the energy density as: 
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where we write 6 (p) as 9 in Eq. |4] and the following Eq. 
[5] for compactness. 

Substracting the energy density of ferromagnetic state 
and multiplying the measure in space integral, we obtain 
the reduced one dimensional energy density ^: 



Anp 



J{P) (\d9_ 
J (00) \2dp 

-13 [cos (0) - 1 

which satisfies dp^ 



su\{9Y K sin (26*) 
4p2 



2p 



(5) 



J tot ^ ^fej- 

To get the optimized function 9 (p) that minimizes the 
total energy, we construct the "imaginary time evolution": 

89 (p, t) _ { 8,^ 8 8^ 
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89 [p) 8p8[8p9{p)] 
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which is a second order partial differential equation for 
function 9 (p) that should be supplemented with appro- 
priate initial condition. The boundary condition is cho- 



9{Q,t) 

0' (Rcnt^r) 



0, 



which means the spins point downward at the center of 
a skyrmion and have the same 0-angle far away from 
the center. The number i?cut is the numerical cut of 
the theoretically infinite system. The initial condition is 
chosen randomly as: 



'(P,0) 



Ttt p = 

I random (0, n) p ^ 0, 



where random (0, tt) generates a random number from 
to TT. We choose i?cut — 10 and slice the space inter- 
val [0, i?cut] to 1000 pieces, then we approximate the 
space derivative by finite difference and do time integral 
using Runge-Kutta method. We numerically integrate 
this initial-boundary-value problem until 9 (p, r) does not 
change with time. The results for different random initial 
conditions are the same, so we find the global minimum 
of the energy functional and the globally optimized func- 
tion 9 (p). 



B. Numerical results 

Now we briefly introduce the phase diagram of a two 
dimensional chiral magnet in zero temperature. For a 
uniform system with constant J and D in space, as we 
increase the magnetic field from zero, the ground state 
changes from helical state to skyrmion lattice state and 
then to field polarized ferromagnetic state [H O . 
In the phase boundary of skyrmion lattice phase and field 
polarized ferromagnetic phase, an isolated skyrmion has 
averaged energy density that is the same as ferromag- 
netic state, so statistically, a dilute population of isolated 
skyrmions is allowed. This so-called skyrmion-gas state 
is the state of interest in the following discussion. It 
is worth noticing that there is a strong resemblance be- 
tween this skyrmion gas and the vortices in the type-II 
superconductor at the lower critical magnetic field H^i . 

Below we show the p-dependence of coupling J, opti- 
mized trial function 9 (p) and energy density ^ in dif- 
ferent parameters (Fig. [T|). We choose n — 1, (3 — 1.6, 
and ^ = 2 for three pinning strength A = 0, A = 0.5 and 
A = 1. We choose ^ = 2 so that the size of the pinning 
center has the same order of magnitude as the radius of 
a skyrmion (so we can make patterned pinning centers 
as dense as possible, see Section. IIICl. It is shown that 
for the uniform case A = 0, the energy density profile 



has a positive maximum and a negative minimum in the 

core region of a skyrmion. If we tune the magnetic field 
/3, there is a /3c ~ 1-6 in which the positive "energy cost" 
and the negative "energy gain" cancel each other (Fig. |2| . 
This is precisely the aforementioned phase boundary of 
skyrmion lattice phase and field polarized ferromagnetic 
phase. 

If there is a local maximum of coupling J, i.e., A ^ 0, 
the core energy is lowered. For the same /3, a larger 
A makes the energy even lower (Fig. [2]), thus there 
can be a negative energy skyrmion in the skyrmion 
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Figure 2: (Color online) Single skyrmion energy versus mag- 
netic field for different strength of pinning center A. It 
is shown that a larger strength of pinning A gives a lower 
skyrmion energy J dp^ for the same magnetic field. So for a 
pinned skyrmion, it is harder to polarize it by magnetic field. 



the following Hamiltonian [TWilET] : 



r.S 



(7) 



Figure 1: (Color online) p-dependence of exchange strength 
J (p) (a), optimized trial function 6 (p) (b) and energy den- 
sity (p) (c) for different pinning strength A. It is clearly 
shown that a local maximum of J lowers the core energy of a 
skyrmion. So a skyrmion is more likely to stay at this pinning 
center. The dashed line in (b) indicates 6 = n/2. 



gas phase (even in field polarized ferromagnetic phase). 
This pinned skyrmion becomes more stable than a "free 
skyrmion" from an energy point of view. This is what we 
call the pinning effect of a skyrmion by a local maximum 
of magnetic exchange strength. Besides lowering the en- 
ergy, the local pinning center also makes the radius of a 
skyrmion a little larger. 



where Sr with — 1 are local magnetic moments at 
positions r = {x,y), x,y ^ - {N - 1) /2, + (A^ - 1) /2 
for even number N, and S — e^jSy are two bond vectors 
in two dimensions. The dynamics of this magnetic system 
subjected to an applied electric current is obtained by 
Landau-Lifshitz-Gilbert (LLG) equation with a current 
term: 



xSr), (8) 



dSr „ , . „ dSr 

— = 5, X H^g (r) -aSrX — 



where is the (5-component of the applied electric cur- 
rent and a the phenomenological Gilbert damping pa- 
rameter. The effective magnetic field -ffgff (r) is obtained 
by the Hamiltonian as: 



III. DYNAMICAL ANALYSIS 

A. Landau-Lifshitz-Gilbert equation 

To study the dynamics of the pinning process, we use 
discrete model of chiral magnets on a square lattice with 
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Figure 3: (Color online) Snapshots of the Landau-Lifshitz- 
Gilbert dynamics of a skyrmion around a pinning center on a 
40 X 40 lattice. The pinning center is at the origin point (0, 0). 
Initially we put a skyrmion at (16, 0), then the offset skyrmion 
undergoes a circular motion around the pinning center and 
finally stays at the pinning center. 



The current term in the LLG equation — ^gjgSr x 
V+(5 ~ '^r-s) ^ '^'^ ^ discrete version of the 
continuous term — (j • V) n |19 |[24H27] . Here we take the 
cross product operation two times to the finite difference 

term (S 



r+S 



r-sj /2 to remove the parallel compo- 
nent with respect to Sr, so that the local spin module 
conserves during the time evoluton. 

We simulate the dynamical pinning and depinning pro- 
cesses as follows: we use an TV by iV square lattice 
with periodic boundary condition, the site-dependent ex- 
change strength and DM interaction 



-^r-S = Jq < 1 + A exp 



(.T + l/2f + (y + l/2f 



r.o 



y2tan(27r/d) 



(10) 



are chosen to have a Gaussian profile, where Jq is the 
discrete version of J (oo) in Eq. [S] and d is the wave- 
length of the chiral magnet on the lattice. Here we have 
shifted the position r = [x. y) in the definition of J to 
the center of the elemental plaquette in square lattice 
(a; -I- 1/2, y + 1/2) so that the periodic boundary con- 
dition is satisfied. We numerically integrate the LLG 
equation using fourth-order Runge-Kutta method and 
follow the position of the skyrmion. We choose time-step 
dt = 0.001 which is small enough for our study and we 
simulate the time evolution from t = up to i = 1000. 
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Figure 4: (Color online) (a) System energy versus time for 
the pinning process. The energy decreases during the whole 
process and finally to a stationary value, (b) The pinning 
path of the skyrmion. The skyrmion moves to the pinning 
center circularly. 



B. The pinning process 

We first choose iV = 40, Jq = 4, A = 0.5, ^ = 8, = 18, 
B = 0.9 e^, and a = 0.2 with zero current. For the initial 
condition, we put a single skyrmion at position x — 16 
and y — which has an offset from the pinning center 
X — y ^ 0. The result of the time evolution is that the 
skyrmion undergoes a circular motion around the pinning 
center with the distance to the pinning center be smaller 
and smaller and finally stays at the pinning center (Fig. 

We have the following heuristic explanation of this pro- 
cess. Because the pinned skyrmion has lower energy than 
a free skyrmion and LLG equation has a damping term, 
the offset skyrmion has a tendency to move to the pin- 
ning center. If we treat a skyrmion as a point particle, 
the LLG dynamics of this single skyrmion has an effective 
kinetic term that gives transverse motion related to the 
pinning force f26|, that is, skyrmion motion is not parallel 
to the "attractive force" of the pinning center but makes 
an angle to the direction of the force. Combining the 
above two factors we get a circular motion with smaller 
and smaller distance to the pinning center. 

In Fig. [4] we show the evolution of the system energy 
(minus the energy of ferromagnetic state) versus time and 
the path of the skyrmion in the pinning process. We can 
see the system energy decreases with time, this is caused 
by the damping term in the LLG equation. 



C. The depinning process 

Now we study the depinning process of a pinned 
skyrmion. Because a pinned skyrmion has lower energy 
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Figure 5: (Color online) Snapshots of the Landau-Lifshitz- 
Gilbert dynamics with current for the depinning process of 
a single skyrmion. We apply current j — (0.2,-0.1) to the 
pinned skyrmion and the skyrmion moves and leaves the unit 
cell from the right boundary and enters the next unit cell from 
the left boundary. Then the current is turned off (at = 200) 
and the skyrmion is attracted by the new pinning center and 
finally is pinned. After the whole process the system state 
returns to the initial one. But in fact the skyrmion is pinned 
by another pinning center because we used periodic boundary 
condition. 
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Figure 6: (Color online) The evolution of the system en- 
ergy versus time (a) and the depinning-pinning path of the 
skyrmion (b). The current is turned on at t = and the 
work done by the current enhances the system energy. Af- 
ter the current is turned off a,t t = 200, the system energy 
starts to decrease because of the damping, and finally reaches 
to a stationary value which is the initial value. After the 
whole process the system state returns to the initial state. In 
this process the skyrmion travels through the boundary of the 
simulated unit cell (red dashed line). 



than a free one, to depin a skyrmion, we need to do pos- 
itive work to it. In order to do work to a spin system, 
we can use a time-varying magnetic field or an applied 
adiabatic electric current. The magnetic field method is 
often used to sustain certain eigenmodes of the spin sys- 
tem near the ground state. According to Mochizuki |14| . 
a time- varying in/out-of plane magnetic field with cer- 
tain frequency can trigger the rotating/breathing inter- 
nal mode of a single skyrmion. In this paper, we choose 
the electric current method to push the pinned skyrmion. 

To simulate the depinning process, we apply a square- 
current-pulse of the form 

j{t)=j,e{t)e{^-t) 



to the pinned skyrmion configuration, where is the time 
when we stop the current and O (t) is the step-function. 

We have adopted periodic boundary condition in our 
simulation of the LLG equation, so we are implicitly sim- 
ulating an infinite pinning center lattice in which each 
pinning center lies at the origin point of our simulated 
N hy N unit cell. We expect the following phenomenon 
for the depinning process. When the current is turned 
on, the pinned skyrmion is pushed by the current to 
move. The direction of the skyrmion motion is deter- 
mined by: (1) the pinning force, (2) the pushing force 
(spin torque) of the current and (3) the tendency of the 
skyrmion to move not parallel to the forces. Let the di- 
rection of skyrmion motion to be in the x-direction. At 
some time the skyrmion would leave the unit cell from 
the right boundary then enter the next unit cell from the 
left boundary and then move to the pinning center in this 
unit cell. Roughly at this time we turn off the current so 
the skyrmion is attracted by the pinning center in this 
unit cell and does circular motion just as the above pin- 
ning process and finally stays at the new pinning center. 

To simulate this process, we use the parameters Jq = 
(0.2,-0.1), = 200 and other parameters that are the 
same as the pinning simulation. We choose the current 
to be not parallel to the a:-direction so that the skyrmion 
moves approximately in the x-direction. The numerical 
result is in agreement with our expectation (Fig. |5|. 

In Fig. ^a.) we show the evolution of system energy 
versus time. The system energy first increases and then 
decreases when the current is applied but always de- 
creases when the current is turned off. This is because 
the current could do positive or negative work to the 
skyrmion depending on whether the current is pushing 
the skyrmion away from or toward a pinning center. But 
it is the positive work done by the current that depins the 
skyrmion. The system energy must decrease without the 
current due to the damping. The path of skyrmion mo- 
tion in the whole process is shown in Fig. [oj^b), where the 
red dashed line indicates the boundary of two successive 
unit cells of the simulated periodic lattice. 
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D. Critical current of the depinning process 
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Figure 7: (Color online) Energy evolution (a) and skyrmion 
patli (b) in different applied currents. For low currents up 
to fi = 0.53 the skyrmion cannot escape from the pinning 
center at the origin point (0, 0) and the system energy cannot 
exceed the barrier indicated by the dashed line at roughly 
H — _ff£g = 4.7. When the current is as large as ^ = 0.54, the 
skyrmion escapes from the pinning center and the depinning- 
pinning process happens several times because we do not turn 
off current in this simulation, so the system energy has a large 
oscillation, (parameters; A = 0.5 and ^ = 8) 
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We found in our simulation that the pinned skyrmion 
would not travel to other pinning centers but instead cir- 
culates around its own pinning position and finally stays 
at a nearby position if the applied current is not large 
enough. So there is a threshold critical current to depin 
a skyrmion. 

To study the critical current of depinning, we do a 
series of simulations with increasing current intensity 
j — /i (0.2,— 0.1) with ^ = 0.51, 0.55 for a pinning 
strength A = 0.5. The result shows that there exists a 
critical current below which the skyrmion cannnot es- 
cape from the pinning center but circulates around the 
pinning center and the energy will not exceed the barrier 
energy set by the pinning center; but when the current 
is larger than the critical current, the skyrmion could es- 
cape from the pinning center and the energy exceeds the 
barrier energy (Fig. [7|. 

Next we simulate the critical current in different pin- 
ning stength A and we find an almost linear dependence 
of the critical current = fi (0.2, —0.1) with respect to 
the pinning strength (Fig. |8|. The approximate relation 
is /i w A. 



E. Estimation of the critical current 

To estimate the critical current, we rescale the LLG 
in real material to our dimensionless form. Let the let- 
ters with a prime represent physical parameters in real 
material. First we do space rescaling, supposing the real 
material has magnetic exchange J' and crystal lattice 
spacing a', then because the magnetic exchange scales as 



the inverse square of length, we have a J 



'-J'. Here 
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J is the rescaled magnetic exchange for our simulated 
space unit defined by 18a « A', where A' is the helical 
wavelength in real material. Next we do time rescal- 
ing, in our simulation we have implicitly chosen the time 
unit to be to = Ah/ J. After these rescalings, the unit 
of the electric current is jo = (Si/p) [qe/ (a^io)] where S 
is the spin quantum number of the local magnetic mo- 
ment, p is the fraction of spin polarization of the electric 
current. Next we substitute real parameters chosen as 
a' = 0.4nm, J' « 3meV, A' sa 60nm, p ss 0.1 and S* w 1 
PT) . then we get jo w 2.4 x lO^A/m . With the above 
result we have jc ~ 5.3A x 10*A/m^ (here A refers to the 
pinning strength with order of magnitude 10^^ in our 
simulation) and the typical pinning and depinning time 
around lOOOio ~ 10ns. 



Figure 8: (Color online) Critical current versus pinning 
strength. The critical current — /i(0.2, — 0.1) needed to 
depin a skyrmion varies almost linearly with respect to the 
pinning strength A. 



F. Multiple-skyrmion depinning 

If there are more than one skyrmions around a pin- 
ning center, each skyrmion will feel simultaneously the 
attractive force of the pinning center and the repulsive 
force of other skyrmions because skyrmions repel each 
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Figure 9: (Color online) Snapshots of two skyrmion depin- 
ning. The current is chosen to be smaller than the critical 
current. But one of the two skyrmions is depinned. This is 
because the depinned skyrmion is subjected to both the pin- 
ning force and the repulsive force of the other skyrmion and 
the net result is the weakening of the pinning effect. This 
figure clearly shows that the depinned skyrmion is smaller 
than the pinned one, which is a conclusion of our stationary 
analysis. 



other. The net result is that the pinning effect is weak- 
ened. We simulate a two skyrmion case with A = 0.5, 
j = 0.4 (0.2, —0.1) which is below the critical current for 
a single skyrmion. The result is that one of the skyrmions 
is depinned by this current (Fig. |9|. Generalizing to more 
skyrmion case, we may say that multiple skyrmion con- 
figuration is more easily depinned. This is a well-known 
experimental result |16| . 



G. Influence of the non-adiabatic spin-torque term 

Besides the current term in our LLG equation, in 
conducting chiral magnets, there is another current re- 
lated term which is the "non-adiabatic spin-torque term" 
[241 E51 [27] . Combining this term to the existing current 
term we have: 



(j ■V)S^{j-\/)S + pSx [{j ■ V) S] 



(11) 



where /3 represents the strength of the non-adiabatic 
torque. It is known that this term alters the trajectory 
of single skyrmions. When there is no pinning center, 
the skyrmion would move parallel to the applied current 
if the value /3 equals to the Gilbert damping a, other- 
wise the skyrmion motion would have a transverse com- 
ponent to the right/left of the current direction if /? is 
larger/smaller than a [27j . 
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Figure 10: (Color online) Skyrmion trajectory for the depin- 
ning process in different /3-values. It shows that the main 
influence of the /3-term is the direction of the transverse com- 
ponent of motion with respect to the current direction. Be- 
cause there exists a pinning center in our simulation, even if 
/3 = Q, the skyrmion still moves in a curved path. 



When there is a pinning center, the /3-term has similar 
effect. In the following simulations we still use the current 
j = (0.2,-0.1). In Fig. 10 we show that the trajectory 




of skyrmion in the depinning process is altered by the 
/3-term. Because there is a pinning center, even if /3 = a 
the skyrmion still moves in a curved path. 



IV. SUMMARY AND OUTLOOK 

We proposed a mechanism to pin skyrmions by a local 
maximum of magnetic exchange strength J in the con- 
ventional Heisenberg model with DM interaction. Physi- 
cally, this local maximum of exchange strength is caused 
by a local maximum of the density of itinerate electrons in 
the material. We studied the effect of this local maximum 
by stationary analysis of the energy functional and found 
that this pinning center makes an isolated skyrmion more 
stable by lowering the core energy of it. Besides the sta- 
tionary analysis, we also studied the dynamics of the chi- 
ral magnet with pinning center by LLG equation with 
a current term. We found that an offset skyrmion near 
the pinning center is attracted to the pinning center and 
undergoes a circular motion around the pinning center 
and finally stays at there. This pinned skyrmion can 
further be depinned by an applied current pulse, which 
pushes the skyrmion away from the pinning center and 
let the skyrmion be attracted by other pinning centers. 
Our simulation shows that there exists a critical current 
below which the pinned skyrmion cannot be depinned. 
We did a series of simulations and found that the critical 
current depends approximately linearly to the pinning 
strength and we estimated the critical current to have 
order of mag nitude 10^ ^ lO^A/m . With the pinning 
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and depinning processes, we actually found a way to ma- 
nipulate the skyrmions in artificially made pinning center 
lattice. This pinning center lattice is expected to be re- 
alized by putting metal grains on the surface of a chiral 
magnetic thin film in a certain way. 

Our pinning and depinning mechanism resembles in 
spirit to the mechanism of the "magnetic domain-wall 
racetrack memory" proposed by Parkin et. al. |28] 
in which they use magnetic domain walls to represent 
information but here we use single skyrmions. The 
main components of a racetrack memory are: (1) a 
writing unit, (2) a racetrack and (3) a reading unit. For 
skyrmions, the writing unit could be realized by the 
method of "skyrmion generation by a circular current" 
[27] . Our work can be thought of as a model of the two 
dimensional racetrack for single skyrmions. In Parkin's 
work they use patterned notches along the edges of the 
racetrack to make pinning centers for the domain walls 
and they use electric current to depin the domain walls 
to make the information "flow". In our present proposal 
we use patterned local maximums of exchange strength 



to pin skyrmions and still electric current to depin 
them. The critical current to depin skyrmions from 
disorder is roughly lO^A/m^ [TB], which is much smaller 
than our estimated critical current to depin skjrmions 
from artificially made pinning center lO^A/m . The 
above two values are both much smaller than the critical 
current in racetrack memory which is of order lO^^A/m^. 
Skyrmions in conducting chiral magnets can be more 
easily depinned and causes much smaller Joule-heating 
for potential application. 
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